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Abstract 

In this paper we intend to extend some ideas of a recently proposed Lorentz-invariant Bohmian 
model, obeying Klein-Gordon equations, but considering particles with a spin different than zero. 
First we build a Bohmian model for a single particle, closely related to the Bohm-Dirac model, but 
using Nikolic's ideas on a space-time probability density and the introduction of a synchronization 
parameter a for the particle trajectory. We prove that, in the case of a spinor whose component 
phases become equal when h 0, '\t is possible to use the usual Bohmian Mechanics procedures 
to obtain the classical limit. We extend these ideas for a two-particle system, based on two non- 
local but covariant Dirac-like equations. 

KEYWORDS: Bohmian Mechanics, Dirac Equation, Nikolic Time, Born-Nikolic Distribution, Non- 
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1. Introduction 

From the Bohmian point of view there is no need for big changes to include, on the theory, 
particles with spin different than zero. In spite of some conceptions about spin, including 
intuitions by David Bohm [1], in further developments of the Bohmian Mechanics it is not 
considered that particles are intrinsically spinning (see Reference [2] for a very good introduction 
to Bohmian Mechanics in a general way including spin). The usual analogy with a little planet 
spinning around its own axis has been fully abandoned by the modern Bohmian theory. 

Bell says that the only hidden variables to be considered are the particle positions. He states 
explicitly that this is a great merit on Bohm's theory. His reasoning is that every measure on a 
laboratory can be reduced to positions. For example, the position of the needle on a measurement 
apparatus or the electron positions when hitting an oscilloscope screen. It is somewhat odd, 
according to Bell, to consider the position of the particles as a hidden variable. Precisely, it is the 
position what immediately reveals itself on a Quantum Mechanics experiment. From Bell's point 
of view, position would be the measured variable and the wave function would be "the Revealed 
Variable". The wave character shows up after doing experiments with a big number of particles 
[3]. 

However, the term "hidden variable" is important in the history of Quantum Mechanics. It is 
usually considered under that term anything being added to the wave function with the purpose. 



for example, of building an Objective Quantum Interpretation, or also, though is not the same, 
recover the deterministic character of Physics on this level. It is from that point of view and only 
from that point of view, that the Bohm Theory is a Hidden Variable Theory. The spin, although this 
term may sound anachronistic [4,5], would be a property of the wave function and not one of the 
particle. Therefore, it cannot be considered another hidden variable. In a Stern-Gerlach 
experiment where you put a magnetic field for the particle to pass through, depending on the 
position wtiere tlie particle is finally detected we say we obtain as a result "up" or "down" [4,5]. 

In Non-Relativistic Bohmian Mechanics a spin-1/2 particle wave function can be thought as a 
function: 

(p: C^, 
the so called spinor. Or, in a different notation. 



(<Pi\ ^ ((piW\ 

\(P2) \(p2(ix))' 



where x is called the generic position of the particle. We will leave the time dependence implicit. 
That is, all the wave functions considered, in general, depend also on time. So on this introduction 
instead of writing (pt as in [4], we will simply write (p. 

Now, for N particles the wave function could be thought as a function: 

On this introduction we will restrict the discussion to two particles for clarity purposes. Then 

(p: m.^ ^ 



or 



<P = 



<P2lix^'\x^'^) <P22ix^'\x^'^) 



Here x^^\x^^^ are the two particles generic coordinates. The (Born) probability distribution for the 
two particles is given by 

<P*ilSijSim(pjm = \(pll\^ + \(P12\^ + \(p2l\^ + \(p22\^ ■ (1-1) 

In the left-hand side term we have summed on the repeated matrix indexes. The particle density 
currents are determined by 

/'^ = ^WilSijSlm^^^^<Pjm - <PilSijSirnV^^^(pJ^} 

(1.2) 

F^=^^[<PtlStjSir,V(^'^(Pjrn-(PaSijSirnV^'^(p;^]. 



On the base of (1.1) and (1.2) we can propose the Bohm-like guide equations for the particles in 
the standard way: 

dt 2imi <PaSijSim(pjm 

dt 2 17712 (P*iAjSim(Pjm 

In these expressions we sum on repeated indexes independently on the numerator and 
denominator. Also on the right-hand side, both expressions are evaluated on the particles 
configurational or actual coordinates ^'^^^ X^^\ Then we prove, in a standard way that if the 
particles are distributed according to (1.1) at time t=0, once this expression has been normalized, 
it will remain distributed that way for every time t. These are very well known results of Non- 
Relativistic Bohmian Mechanics. 

Let us suppose, as Bell does [3] the particular case when q)ii is separable that is cpy (x^^\x^^^) = 
Xi(f(^))f/(x®).Then 

dt 2tmi XHl^ij^lrnXj^m 

h (^!Si^^m)x!StjV^'hj-i^lSlm^m)XiSijV^'h; 

2^mi ix*StjXjMSirn^m) 
Now ^tSim^m = ^Am^m since it is real. Therefore 

djW^ h xtSjjV^'hj-XiSijV^'hj ^ h Imjx^V^'^x) 
dt 2imi XtSijXj rui x^X 

and analogously 

dt ■ 

This shows the difference between a separable wave function and a not-separable one in Non- 
Relativistic Bohmian Mechanics, following Bell's ideas [3]. 

In Sect. 2 of this paper we extend this kind of ideas to Special Relativity. We base the study on the 
Dirac equation, the relativistic equation for the electron. We use Minkowski notation for Special 
Relativity in which it is not needed to distinguish between covariant and contravariant indexes [6]. 
On this notation the coordinates of an event are denoted by x, ict. In this case the wave function 
is a four component spinor. We also assume that the particle is interacting with an external 
electromagnetic field, measured in Gaussian units. As it is usual in Bohmian literature, we 
supplement the wave function with a Bohm-like equation that establish in a precise way how the 



particle moves. We want the Bohmian velocities and the whole model to be covariant. We then 
make a proposal for a guide equation considering the following ideas by Nikolic [7]: In the first 
place we consider an space-time probability density instead of an spatial density. We introduce an 
equivariance (quantum equilibrium) condition by introducing a scalar parameter a with time units 
[7,8]. We assume that the probability distribution at time c = is given by our space-time density. 
From this it is proven that the particle probability density will remain so for every c. Since the 
wave function, the spinor, does not depend on a, we have a stationary flow with respect to this 
parameter with a space-time distribution that is also a scalar [7]. 

The only way to prove, for now, that the model presented here is physically reasonable is to study 
the classical limit [9]. We do this on Sect. 3. We do not consider the classical limit in general, but 
just in a case we think it is particularly important. Under a certain hypothesis on spinor, that its 
components phases being equal on the classical limit. We can prove on this hypothesis that: {a} 
The interval da tends to the proper time of the particle dx. This is analogous to the study of the 
classical limit in [10] for a A/-particle model based on the Klein-Gordon equation [7,10]. (h) The 
classical movement law of the particle obtained is the one for classical Einsteinian Mechanics for 
one particle subject to the Lorentz Force, (c) Two kinds of classical solutions are obtained in 
different regions of space-time. One with a negative charge e, the electron charge, the other one 
with a positive charge -e. 

We believe that the hypothesis of phase equality on the classical limit is not entirely arbitrary. This 
hypothesis tells us that the resulting classical system has a well defined action. Let us suppose a 
fixed action 5 (function of the space-time coordinates) that fits as the classical limit of the 
component phases of a set of spinors. It can be proved that this set is a subspace of the solution 
space of spinors from Dirac equation. We conjecture that the Bohmian problem based on Dirac 
equation, restricted to this subspace, is particularly physically relevant. We believe that a very 
similar hypothesis can be used for the classical limit following the method of [9] on the Non- 
Relativistic version of our model, Pauli equation from a Bohmian point of view. 

In Sect. 4 we study the model extension for two particles. The particles interact with an external 
electromagnetic field but there is no electromagnetic interaction with each other. It is possible to 
build a model for N particles (Some kind of ideal gas for particles with spin different than zero [7, 
10]). For the sake of simplicity we will work with only two particles. We use the Nikolic's idea [7, 
10] of a density on the space-time configuration space of the two particles x^^^, ict^^^ , x^^\ 
ict^^\ As it can be seen the model is multitemporal. The times ict^^\ ict^^^ are treated on an 
equal footing with spatial generic variables. The Wave Function, a 16-component spinor, is 
supplemented with Bohm-like equations. These equations determine how the configurational 
variables X^^\icT^^^ , X^^\ icT^^^ evolve in terms of a parameters. This parameter, with time 
units, is a scalar and represents an extension of [7] to a model built on Dirac equation. The model 
is Non-Local and at the same time Covariant and some ideas of Bell [3] are analyzed on the Special 
Relativity context. 



Finally, in Sect. 5, we present our conclusions. 



2. A Modification of The Bohm-Dirac Model. 



In this paper we will develop a Bohm-like theory for the Dirac wave equation in the one particle 
case. We start with the Dirac equation: 

Y^{p^-^\)^ = -mcxp (2.1) 

where 




the spinor, plays the role of the wave function for one particle [11-12]. On the other hand, 
represents an external electromagnetic field, evaluated on the particle generic coordinates, e is a 
negative charge equal to that of the electron, and the matrixes are defined in terms of the 2x2 
Pauli matrixes as follows: 
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l^ ro -i 
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Usually it is defined j/i = j/;"*")? where j/;"*" = (j/;i* ip2* V's* '/'4*)- 
In this case 

is real [11-12]. The expression (our candidate for space-time probability density) 

\M (2.2) 

is a scalar. We use absolute value because we want a positive definite expression for the space- 
time probability density. The term ijjilj could easily be negative on certain regions of space-time. 
In the first places/; = Sif), where 5 is the transformation matrix that transforms the spinor from 
one frame of reference to another. Therefore ip'^ = xp'^S'^ and ip'^P = xp'^S'^p = ■[p^fiS'^. Here 
we have used the identity 5"*" )S = ps~^. By multiplying both sides of the identity b'^xp = Sxp we 



have = j/'"'"/?S ^Sip = Or = and considering the absolute value on both 

sides we obtain that \4>'4''\ = This tells us (2.2) is in fact a scalar. 

The following are very well known facts. If t/; satisfies the Dirac equation then 

cxpYn^P (2.3) 

is a four-vector and 

^ (cV;r^V) = 0. (2.4) 



From a Bohmian point of view Equation (2.1) does not completely describe the dynamics of the 
particle. We have to supplement (2.1) with a Bohm-like equation that describes how the velocity 
of the particle is determined by the spinor xjj. We propose 

where the right-hand side of the equation is evaluated on the particle configurational coordinates 
X,icT. On the other hand, da is a parameter with time units that we will use to describe the 
dynamics (Let us call da "Nikolic time for the model", since we are using a very similar idea to the 
one used by Nikolic in [7] to build a Lorentz-lnvariant Bohmian Model, based on the Klein-Gordon 
equation), da is clearly a scalar since and cipYii4'/\^^\ are four-vectors. Let's rewrite (2.4) as 

V,) = 0. (2.6) 



On the same way than in Non-Relativistic Bohmian Mechanics, the system dynamics must be 
supplemented with a probabilistic hypothesis on initial conditions. Let us assume that the 
probability that the particle position four-vector is in the space-time volume AxAyAzcAt 
surrounding on cr = is proportional to p(x^) AxAyAzcAt = |i/^i/^|AxAyAzcAt. Thus, as a 
consequence of the equations (2.5) and (2.6), the particle position four-vector will continue to be 
distributed according to the same probability density [7, 10]. Since we are using a very similar idea 
to Nikolic's [7] in the construction of the model, we will call \ipip\, once normalized, Born-Nikolic 
probability distribution for the model. 

In Special Relativity the volume of a region of space-time is independent on the frame of reference 
used for its calculation. This condition is given by 

d(x' ,y' ,z' ,ct') _ ^ 
d{x,y,z,ct) 



It follows that the volume of a region of space-time is also Lorentz-lnvariant [10]. The probability 
assigned by the Born-Nikolic rule (for some value of a) to a given region V of the space-time for 



the particle is given by the integral of p(x^) = \-^rp\ on V. It follows that the probability assigned 
to a certain region V in an inertial frame of reference will be the same than the probability 
assigned to the Lorentz-transformed region V in a different frame of reference. In other words, 
the calculation of the probability is independent on the frame of reference used. Therefore, if the 
particle position four-vector is distributed according to the Born-Nikolic rule proposed here in an 
inertial frame of reference for some value of a they will be so distributed for all inertial frames of 
reference and for all values of a. The "quantum equilibrium condition" is thus, in the model, 
Lorentz-invariant [7, 10]. 

A very important hypothesis that we are considering, not so much in the construction of the 
model, but in the probability considerations is that p(x^) = is normalisable. This is not 
proved in any way. Nevertheless, let us consider for a moment, as an exercise, xp'^ilj as a space- 
time probability density (This disagrees with Nikolic's proposal [7, 10] because this quantity is not 
scalar but the fourth component of a four-vector). Now, the integral ofip'^ip with respect to space 
is time-independent due to the Dirac continuity equation (2.4). If after that we integrate respect to 
time the integral diverges and so ifj'^ifj is not normalisable if we think of it as a space-time density. 
However this argument cannot be applied to \ipip\ because the spatial integral of this quantity is 
not constant on time. We will consider the provisional hypothesis that there is no problem with 
normalization of this function in space-time. 

In the following we consider the relations fSyi = i = 1,2,3 and = i- It is important to 
compare the exposed ideas and some important ideas of Bohm's school about construction of a 
Bohmian Model for Dirac equation. Our equation (2.5) can be decomposed on spatial and 
temporal parts as follows 

It follows that dT Ida = which implies 

dX cxfj'^dilj 



dT xp+xp 



(2.7) 



This constitutes the Bohm-Dirac theory as formulated by David Bohm [13-14]. We can consider the 
right-hand side of (2.7) as a function of the configurational coordinates X and timeT. Then the 
wave equation can be written as 

dip ^ r ^ e ^\ , 
ift— = ca- [p — - Aj xp + (ecp + m.c P)xp. 

Here the spinor is seen as a function of the generic coordinates x of the particle and 
configurational time T. The continuity equation can be written as 



— (ci/^+a,!/.) + — (1/^+1/.) = 0. (2.8) 

Then if at time 7 = the particle is distributed with probability V^V'/ once normalized the 
function in the space, this will remain so for every T. This follows from the continuity equation 
(2.8) and the Bohm-Dirac equation (2.7). This property of "quantum equilibrium" is only valid on a 
preferential frame of reference. That is, there's quantum equilibrium in a preferential frame of 
reference but not in any frame of reference. 

Since what we are interested, like Nikolic [7, 10], in a space-time density in the relativistic context 
we will work with the guide equation (2.5) and the space-time density rather than with the 
guide equation (2.7) and the space density i/^"*"!/^. Our idea is that this procedure (insisting on the 
space-time probability density) allows us to extend the model to a Non-Local model for N particles 
that remains Lorentz-lnvariant. 

3. Classical Limit 

We now intend to extend the ideas of [9] to obtain the classical limit. We will not do it in general 
but only for spinors that satisfy the following hypothesis. 

Hypothesis of phases equality on the classical limit. Let xpi = R^expiiSi/K) be the i-th component 
of the spinor. Let us assume that 

lim Si = 5 

where S is independent or\i {i = 1, 2, 3, 4). In this section and only in this section we will use S for 
that limit. On the limit when h tends to the phases of the different components of the spinor 
become equal to each other. If we write 5j as a series in h we obtain 

Si = s + sl;^h + s^^h^ + sph^ + ■■■ 

where S is independent on the index i. Observation: = [h], that is, every term on the 

series has units of action. Always that the spinor has the form just described we will say that we 
are using the hypothesis of phases equality. Let us consider the set of spinors with the property 
that the hypothesis of phases equality is satisfied on the limit when ft -» in such a way that 5 is a 
fixed function. It can be proved that this set is a subspace of the space of all spinors that satisfy the 
Dirac equation. We may consider a Bohmian dynamics restricted to this subspace of spinors. As a 
clarification, if we add two spinors in different subspaces of this kind the sum spinor does not 
satisfy the hypothesis of phases equality on the classical limit. 

_.s 

For calculation convenience, let us define \p'^ = e '^i/^. From the hypothesis of phases equality on 
the classical limit we can define the "classical spinor" as 



_.5 

^class _jjj^^c_ ijj^g '^hip. 
ft-»0 ft-»0 

The "classical spinor" is only correctly defined when using a hypothesis as pointed out before. This 
spinor does not depend on h and its components are ijjf'^^^ = R^exp {iS^^^^. 

Now, if i/j satisfies the Dirac equation then satisfies the equation 

[Vn - -^ii)\Pii - -^11)^ + -^^yiiYvPiiv ^ = -m^c^^- (3.1) 

The reciprocal is not true. Equation (3.1) has more solutions than the Dirac equation [11]. 
Equation (3.1) is very similar to the Klein-Gordon equation except for an additional term that is in 
a certain way the relativistic extension of the Stern-Gerlach term of the Pauli equation. Let us 
rewrite Equation (3.1) making explicit every operator: 

, d^ip ehdA^ 2eh dtp ieh ^ , 

Let us consider the following facts, always under the hypothesis of the phases equality on the 
classical limit (Appendix 1): 

^2''' dS dS ^i^class 



i) lim f-fi^e"^-^^^) = — — i/; 

h->0 y dx^dx^j dx^ 

h->0 \ ic dx^ ^ ) 

.... ,. / 2eh -i- . dip\ 2e dS , , 

vi) lim f— m^c^e~'ftj/; 1 = —m^c^ib' 



\class 



,2 ^2,i,class 

— Ill, I, c "u/ I — — (/ 
/ 

Therefore 

dS dS , 2e dS , , , 

,1, class A ,,1, class j a a ,,1, class ,yy,Z ^Z,,i,class 

dx^dx/ cax/"'^ +^2^"^"'^ - mcxiJ 

Given that we suppose 0, after some rewriting we have: 

dS e^\/dS e 
c "JXdXi, 



Which tells us that 5 satisfies the so called Hamilton-Jacobi Relativistic Equation. Thus the first 
step of our argument is complete. 



Let us study the guide equation (2.5). Rewriting the Dirac equation on the following way: 



Yii I 1" ^ — e ^0'' Ai,e tixjj'^ = —mce f^xb'^. 



The limit of the last expression multiplied by e when ft ^ is 



\mdXn mc 



A^, )ip 



.class _ 



-cxp 



class 



(3.3) 



Let's use the following notation U,, = —-^ —Ai- 

I* m oXfi mc 1^ 



Then (3.3) can be written as a matrix equation: 



iU^ + c Us Ui - il 

iU4 + c Ui + iU2 -U3 

-U3 -Ui + iU2 -if/4 + c 

^-Ui-iU2 f/3 -if/4 + cy 





/^class \ 
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^i.class 








^i.class 








\,l,class i 







This equation has non-trivial solutions since the determinant is zero. The solution space is a two 
dimension space, so we write the components of the classical spinor in terms of xp^"-^^ , j^^;ass 

^chiss =R^e'ri, ^p^^^ =i?2e ^, ^z'^' ^ j^j ^ 

— iU^ + c 

-iU4 + c 



So we can build the analogous ofxp^ but using the "classical spinor" 



^class ^class _ 



- 2icU4 + c2 



(3.4) 



Since = — ^ — —0) is purely imaginary xjjdass -i^ciass pgg| j^-^^ should be so because 

we proved that is real. Now, the classical version of the current cxpYn'^ that can be calculated 
from the expressions for the "classical spinor" components is 



C-(Jjclass Y^-(Jj 

The expressions (3.4) and (3.5) lead us to 



class = 2icU,iRi + Rj) _ 
U4 + ic 



(3.5) 



Let us evaluate (3.6) on the configurational coordinates^^. We have: 



— = lim — = lim I - I = , , - — -. 

OJ h^O da h^O ^ijjijj^ \^cl(iss ^class \ 

Here dr is the classical limit of da. In other words lim^j^o da = dr. 
We can consider two cases: 
Case (i): 

dXi, _ _1 dS e 



if is in a region where ipdass -fpciass ^ q 
Case (it): 



dXr. 1 dS e 



if is in a region where ipciass -ipciass q 

dA'u 1 dS e . 
Case (0 = U., = A... Let us start writing (3.2) in the following way 

dr m dX,, mc ^ 



1 dS e \fl dS e , 



mdx^i mc ^J\mdXn mc 

If we evaluate (3.9) in the particle configurational coordinates X,, we have ^^E^^ii = or 
equivalently d^^d^^ = —c^dz^. The last expression reveals that dr is precisely the particle 
proper time. Therefore da tends to the particle proper time when ft -» 0. This result is very similar 
to the one obtained in [10] for a Bohmian Model built for a system of particles based on the Klein- 
Gordon equation [7]. If we derivate (3.9) with respect to the generic coordinate we obtain: 

f dS e \ d / dS e 
\dXn c ^ ) dXy\dx^ c ^ 

fdS e ^ \ d^S _e/dS e ^\dA^ 



\^Xf^ c ^ ) dXf^dXy c\dXi^ c ^ ) dx^ 



Then it is proven, following a simple algebraic procedure, that 



If we now evaluate this identity in the particle configurational coordinates Xf^ we obtain that 

d^v d e 



" = U„—U,, = —U„F„, 



dT^ "dX^ " mc " 

d^X^ _ e dX^ 
"^'d^'c'df^""- 

Which constitues the Einsteinian equation for one particle ruled by the Lorentz force in Classical 
Relativistic Mechanics. That is, the ideas described by [9] produce the correct classical limit when 

we assume the hypothesis of the phases equality on the classical limit. The physical meaning of this 
hypothesis is that when we build the classical dynamics resulting from the limit when ^ -» on 
the original Bohmian system, the resulting physical system has a well defined action S. 

Casein) ^=-U, = --— + -A, 

In this case, when we take the classical limit, following a completely analogous procedure we 
obtain 

d'^X^ _ e 

The essential difference is the change on the sign of the particle charge. With this ends the 
derivation of the classical limit for a Bohmian particle associated to the Dirac equation. 

Summarizing, equation (3.2) along with equation (3.7) or (3.8) are the classical counterpart of the 
wave equation (2.1) and Eq. (2.5) (always assuming the hypothesis of the phases equality on the 
classical limit). In both cases we obtain adequate classical limits and the parameter da tends to 
the proper time dr. The latter is entirely analogous to the case where the classical limit is taken on 
a Bohmian Model for the Klein-Gordon equation [7] in [10]. 

4. Two Particles. 

For the sake of simplicity we will consider only two particles. The problem, however, can be 
extended to any finite number of particles N. This would be sort of an ideal gas where particles do 
not interact electromagnetically with each other (for some models in the bohmian literature [IS- 
IS]). In the model, however, the particles interact with an external electromagnetic field. The 
classical limit of a system with more than one particle will not be analyzed in this paper. We 
believe, however, that the hypothesis of the phases equality on the classical limit works for many 
particles as well. 



The wave function for two particles has the following form: 



The Pauli principle in terms of tiiis expression would be 

{4'^4'^) = -Aj {4'\x^'^). (4.1) 

The quantity 4'a*l^iml^ij4'jm obtained by summing on repeated matrix indexes is a scalar and 
equal to 

^iCPlmPij^jm = IV^ll I' + \^12 I' - 1^13 I' " 1^14 P + IV^21 P + 1^22 I' " 1^23 P " IV^24 l' 
- \^3l\' - \^32\^ + \^33\^ + \^3a\^ - " \^42\^ + l'/'43l' + \^A4\^- 

This quantity is real. Every term of the expression is real. Of course this quantity could be negative 
in some domain of the space-time configurational space for the two particles. 

We will try to extend Nikolic's idea. The wave equation would be formed by the following 
equations. In the first one the operator p^^-* acts on the first particle generic coordinates x^^ and 

(2\ (2') 

the operator on the second equation acts on the second particle generic coordinates . 
^4^'-* is the electromagnetic field evaluated on the i-th particle generic coordinates; i = 1,2. 



We believe they are very similar to a model presented in [15] for many particles. On the other 
hand the Bohm-like guide equations would take the following form: 

It can be proved that both of these expressions determine four-vectors. On the right-hand side of 
the expressions (4.2) we sum over the repeated matrix indexes independently on the numerator 
and denominator. The right-hand side of both equations is evaluated on both particle 
configurational coordinates X^^\ icT^^\ X^^\ icT^^\ Here da and \^ii*PimPij4'jm\' once 
normalized, are respectively Nikolic time and Born-Nikolic probability distribution for this model. 



The parameter da is a scalar since everything else in (4.2) is covariant. To prove this (Appendix 2), 
they can be used the following facts: 

(i) ^'sr = ^siS^jipij that shows how the wave function transforms when passing from a 
frame of reference to another and S is the transformation matrix that transforms one 
spinorfrom a frame of reference to another. 

(ii) S'^P = PS~^. This fact was already used on the first part when we proved that is 
a scalar. 

(Hi) Yu = T^^'^Yri^- This relation is very important to prove that the Dirac equation for 

ox II 

one electron has the same form in any inertial frame of reference. 

In [15] very interesting characteristics of this kind of models (Synchronized trajectories) are 
pointed out, even if the author perceives them as a flaw. 

It may be interesting to do a little exercise in order to compare the model with Bell's reasoning [3] 
about one of his models. This reasoning was discussed on this paper introduction for the Non- 

Relativistic case. Let us consider a region of the space-time configuration x^^.ict^^^, xj^\ict^'^^ 

that we will call S where '^Pa* PimPij'^jm > 0. Then if the configurational coordinates ^'^^^ icT^^\ 

are in this region £ the equations (4.2) have the form: 

t* y 

da Tpil*PlmPijTpjm ^^2) 

dxl,^^ _ ciPii*l3ij(l3y^)^^iPj,^ 

These equations are very similar to the ones presented by J.S. Bell in [3] and in general represent a 
non-local model that, however, has the remarkable characteristic of being covariant. Let us 
assume now that the spinor is separable xpn = (piXi (here we are leaving aside the Pauli Exclusion 
Principle (4.1)). Then 

dx!;^'^ C^aPlm{.PYn),,-^im C(Pi*Xl*Plm{pYn),.(Pj*Xm* C(Pi*{pYn),,(Pj* 
f* y *j *j 

d(T ^il* PlmPij^jm (Pi*Xl*PlmPij(PjXm <Pi*Pij(Pj 

and analogously 

d<y Xl*PlmXm 

This kind of argument is used by Bell [3] to prove the difference between Bohmian dynamics when 
the particles are not separable and when they are. They illustrate the concept of non-locality in 
particles with spin. So we believe that with this kind of procedures it is possible in principle to 
extend much of Bell's ideas, about Bohmian dynamics, to Special Relativity. 



The following continuity equations are a direct consequence of Dirac-like equations on which our 
model is based (Appendix 3): 



(1) {c^aPlm{PYn)^j^jm) = 



^ {cral3ijilSY^)^^xlJjm) = 0. 



(4.4) 



Then 

d 



d 



(4.5) 



Let us add equations (4.5) and consider that p (xj^^\xj^^^^ = \ilJiiPimPij'^jm \ does not depend on 
a. We can write an entirely analogous equation to the one presented in [7, 10], that is, a 
relativistic conservation equation for the two particles systenn under consideration: 



(1) f2") 

Therefore if for u = the particle configurational coordinates , are distributed with 
probability distribution p ("equivariance condition") these will remain so distributed for every a. 
Now, it is satisfied the following equality for the jacobians 

d (xW',yW',zW',CtW') d (x(2)',y(2)',z(2)',ct(2)') 

due to the Lorentz transformations [10]. The measure for an element with a fixed volume dV on 
the configurational space, for two particles with generic coordinates xji^\xj^\ is independent on 
the frame of reference used for its calculus. Let us integrate p over an arbitrary volumen of the 
two-particle configuration. Since the element of the 8-dimensional volume is a scalar and p is also 
a scalar, the calculation of the probability ofX^ \X^ ^ being on 7 at a given time a is independent 
on the frame of reference used for the calculation. If the two-particle configurational coordinates 
z; , z; are distributed according to p (once normalized) at time a = they will remain so for 
every a and for every inertial frame of reference. We do not see any reason why these arguments 
could not be extended to a larger number of particles. In that case the spinor would have the form 



where N is the number of particles. 



5. Conclusions 



A Bohmian dynamics associated to the Dirac equation has been built, initially, this has been done 
for one particle and subsequently for two particles. This process can be easily extended to any 
finite number of particles. In the construction of these dynamics Nikolic's idea of a space-time 
probability density (or in general a space-time configuration) has been used and c as a scalar 
parameter for the dynamics has been introduced. In the example for one particle the classical limit 
following the ideas of [9] has been built. This was done only for a certain kind of spinor, which has 
the characteristic that its component phases have the same limit when h ^ 0. We consider that 
this example is not only special but very relevant since an adequate classical limit was obtained. 
This problem needs further study. The example for two particles is interesting because it is non- 
local and covariant (which is a characteristic of the model built by Nikolic [7, 10] for a Bohmian 
Model based on the Klein-Gordon equation). Finally a very similar analysis to the one done by J.S. 
Bell for separable and non-separable Bohmian dynamics is obtained [3]. 
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Appendix 1 

In the following steps we show a calculus to study the limit of (;j and (Hi) (Section 3) under the 
Hipotesis of phases equality in the Classical Limit: 
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Therefore, the following limit is obtained: 
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The limit (Hi) will be considered in the following steps: 
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The limits (ii), (iv), (v) y (vi) can be established in a more direct way. 
Appendix 2 

It will be proved that the particles's currents 

are four-vectors (Section 4). The current for particle 1 in another frame of reference is : 

4^^' =Clpp*ql3qril3Yr,)^p'sr 

To write (2) the following relationship has been used: 

If this expression is conjugated = S*iSgTn^im — ^im^ip^mq 's obtained. This last equality 
has been used to in order to obtain (2). In the following step the expression (2) is rewritten using 
the relation 5+)S = pS''^: 

Jr] ~ ^4'lm^lp^mq Pqr l^pn{Yr]^^^^si^rj4'iJ 
= (^'^Im (-^ip Ppn ) i^mq Pqr ) (x?? )^ ^si^rj ^ij 
= (^^Im {Plp^pn){Pmq ^qr^ijr))^^si^rj ^ij 

The following step is to change the order on the terms in this expression. Then, the relation 

Ya = Y^^'^y^^ ^^^^ (^'^'^ equality is very important in proving that the Dirac equation has the 
same form in every frame of reference) : 

= C^lm Pip Pmq (^n (Kt? )^ -^si ) (^^/^r; ■ (3) 

By multiplying (3) for ^ the following identity is obtained: 

-^JIt^ = (^'^Im Pip Pmq { ^ -^pn {Yt^ ) ) S^j Xpg 



Sumarizing: 



and analogously: 



,(1) ^ ^,(1)' 



ji2) ^ d^r, (2)' 

ax/" ■ 



Therefore, the particles's currents are four-vectors. Next, the expression ^aPimPij'^jm 's 
considered. It'll be proved that is a scalar. 

^pq Pqr Pps^sr ~ ^Im^lp^Tnq Pqr Pps^si^rj^ij 
~ ^Im iplp Pps ) ("^JTiq Pqr ^^si^rj ^ij 
~ ^Imifilp^ps ^^Pmq^qr ^^si^rj 
~ ^ImPlpPmq ("^ps ^si){^^qr ^rj^^ij 
= ^lm Pip Pmq ^pi ^qj ^ij 

Therefore it follows xp'^q PqrPps^sr = ^imPuPmj ^ij and this expression is a scalar. The expression 
\ipimPiiPmj'^ij\i the Candidate to be the probability density distribution in space-time 
configuration for the two-particle system, is also a scalar. Therefore, the bohmian velocities: 

dX^^\ cl|JtmPmjiPY^,)n4'ij 
\Vlml^liPmj^ii\ 

are four vectors. 
Appendix 3 

In the following expression let i to be a fixed value chosen between 1, 2, 3 and 4: 

^Ii{Py,),%. 



This expression is a function of tlie generic variables x^^^ and For a wliile x^^^ will be 
considered fixed. Once I and x^^-* are fixed, the expression j/;^^ will be a column of the spinorthat 
satisfies the Dirac equation: 



Then, it can be concluded that the following continuity equation is true: 

This is a very known consequence of the Dirac equation. Like a consequence of the nature of )S: 
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The following expression can be developed like: 

Each term of the right hand in this identity satisfies the continuity equation (1). It follows 
immediately that: 

-^UtmPmliPY,\lpjl) = 0. (2) 

This is in perfect accord with the definition of the particle current for particle 1. Analogously 

d 



{^lmPij{PY,)^^^jl) = Q- (3) 



These results are in agreement with the definition of particles' currents (Section 4). 



